Fermions in 2D Optical Lattices: Temperature and Entropy Scales for Observing 

Antiferromagnetism and Superfluidity 
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One of the major challenges in realizing antiferromagnetic and superfluid phases in optical lattices 
is the ability to cool fermions. We determine constraints on the entropy for observing these phases 
in two-dimensional Hubbard models. We investigate antiferromagnetic correlations in the repulsive 
model at half filling and superfluidity of s-wave pairs in the attractive case away from half filling using 
determinantal quantum Monte Carlo simulations that are free of the fermion sign problem. We find 
that an entropy per particle ~ In 2 is sufficient to observe the charge gap in the repulsive Hubbard 
model or the pairing pseudogap in the attractive case. Observing antiferromagnetic correlations or 
superfluidity in 2D systems requires a further reduction in entropy by a factor of three or more. In 
contrast to higher dimensions, we find that adiabatic cooling is not useful to achieve the required 
low temperatures. We also show that double occupancy measurements are useful for thermometry 
for temperatures greater than the nearest-neighbor hopping. 

PACS numbers: 71.10.Fd, 37.10.Jk, 71.27.-fa 



An exciting new development in cold atoms is the abil- 
ity to realize in the laboratory simple models of strongly 
correlated fermions in optical lattices P, |3, H, S S @1- 
These studies are motivated by their relevance to spec- 
tacular phenomena in condensed matter physics, like 
high Tc superconductivity, that are not fully under- 
stood. The best known model is the fermion Hubbard 
Hamiltonian[3, Q that captures the physics of antiferro- 
magnetism and, at least qualitatively, d-wave supercon- 
ductivity in two dimensions (2D). The Hubbard model is 
well understood in one dimension (ID), using exact solu- 
tions and bosonization 0] , and also in the limit of large 
dimensions, using dynamical mean field theory (DMFT) 
[lo| . The two dimensional problem, of direct relevance to 
layered high Tc superconductors, is the least well under- 
stood theoretically. New insights into the 2D Hubbard 
model can come from cold atom emulators, given their 
high degree of tunability (interaction strength, chemical 
potential) and absence of disorder and material compli- 
cations. 

The principal challenges for the optical lattice emu- 
lators are to cool down to sufficiently low temperature 
to see interesting phases and to observe the characteris- 
tic order. In this paper we present detailed quantitative 
results from quantum Monte Carlo (QMC) simulations 
^ of (i) the 2D repulsive (U > 0) Hubbard model at 
half filling and (ii) the 2D attractive {U < 0) Hubbard 
model at any filling, that are both of direct relevance to 
ongoing experiments. The reasons for focusing on these 
systems are threefold. Determinantal QMC simulations 
are free of the fermion "sign problem" in both cases 
and we can obtain non-perturbative results on finite sys- 
tems without any approximations. These two systems 
exhibit phenomena of great interest: strong antiferro- 
magnetic correlations and Mott physics for C/ > 0, and a 



Berezinskii-Kosterlitz-Thouless (BKT) transition to an s- 
wave superfluid with a pairing pseudogap above 
Tc for U < 0. Finally, both these models are realizable 
in cold atom systems where, e.g., a Feshbach resonance 
can be used to change the sign of the interaction. 

Our main results are: (1) we determine the charac- 
teristic temperature scales for spin, charge and pairing 
correlations and the superfluid Tc as functions of the in- 
teraction strength. (2) We determine constant entropy 
contours in the temperature-interaction plane, that give 
direct information on how much the system has to be 
cooled, i.e., how low the entropy should be in the ex- 
periments, in order to see interesting physics. (3) We 
comment on the difficulty of using adiabatic cooling in 
2D. (4) Finally we point out the temperature range in 
which a simple observable like double occupancy can be 
used for thermometry. 

The single band Hubbard Hamiltonian is 



(1) 



where da is the fermion destruction operator at site i 
with spin a, ni„ = c\^Cia is the density of fermions with 
spin cr, and Ui — n^o-. For cold atoms a =|, [ labels 
two hyperfine states. We consider near-neighbor hop- 
ping on a square lattice with the kinetic energy e(k) = 
— 2i(cosfca,a -I- coskyo). In the determinantal QMC cal- 
culation, we work in the grand canonical ensemble and 
tune the chemical potential n to obtain the desired den- 
sity p = '^i{ni)/N on an A''-site lattice. The parameters 
t and U can be directly related to the lattice depth 
Vb, tuned by the laser intensity, and to the interatomic 
interaction tuned by a Feshbach resonance. We work in a 
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FIG. 1: Phase diagrams of the 2D Hubbard Models with 
(a) repulsive {U > 0) interactions at half filling and (b) at- 
tractive {U < 0) interactions away from half-filling, showing 
curves of constant entropy (per site) s = S/Nks, obtained 
from QMC simulations on A'' = 10^ lattices. In (a) Ach is the 
scale for the formation of a charge gap, Amf is the mean field 
charge gap (see text), and Tspin is obtained from the peak in 
x{T) in Fig. 2(b). In (b), the pairing pseudogap tempera- 
ture T* is obtained from the peak in xiT) in Fig. 2(c), and 
the Berezinskii-Kosterlitz-Thouless transition Tc is estimated 
from the superfluid density, Fig. [2jd); Tc as a function of U 
also shown in the inset. 

parameter regime where only a single band is populated 
in the optical lattice. All energies and temperatures are 
measured in units of t and distances in units of a. 

The results for the characteristic temperature scales 
(open symbols) and constant entropy curves (filled sym- 
bols) in the (T, C/)-plane are shown in Fig. [1] We first 
describe the results for the repulsive {U > 0) model at 
half-filling p = 1 and then turn to the attractive {U < 0) 
case at arbitrary filling p 1. 

Repulsive Hubbard model at half-filling: Coming 
down from high temperatures the first scale encountered 
is the energy gap to charge excitations. A lower scale 
corresponds to the development of spin correlations. The 
Mermin- Wagner theorem precludes a finite temperature 
phase transition in this 2D system and both these scales 
are crossovers. Nevertheless, these scales (summarized in 
Fig. [1]) have very clear signatures in physical observables 

(Figs. HSU). 

Charge gap: The charge gap is seen in Fig. HJa) where 



we plot the density p as a function of the chemical po- 
tential p.. Despite large error bars, arising from the sign 
problem away from half-filling, we clearly see that the 
compressibility dp/ dp is very small; it should vanish at 
T — 0. The region \p\ < Ach for which dp /dp ~ gives 
an estimate of the charge gap Ach. The t^/t-dependence 
of the charge gap is shown in Fig. [TJa): for T :S> Ach 
we expect the system to look effectively gapless while for 
T <C Ach there is a gap to fermionic excitations. 

For large U the charge gap scales linearly with U and 
is called the Mott gap, but it becomes (exponentially) 
small for low U. To understand this f7-dependence we 
use a simple T = mean field theory (MFT) of spin 
density wave (SDW) ordering that shows the system is 
an insulator at p = 1 for any U > 0. The [/ = 
metal is unstable due to Fermi surface nesting, and the 
ground state is an insulator whose gap is obtained from 
N~^J2i,[4 + ^mf]"^''^ = ^/U- The mean field charge 
gap is Amf ~ texp —2TTy/t/U for U/t <C 1 and smoothly 
crosses over to the Mott gap Amf — U/2 for U/t ^ 1. 
This evolution lj| from a SDW insulator to a Mott insu- 
lator is the U > analog of the BCS to BEC crossover in 
attractive Fermi systems. We see from Fig. [Ua) that the 
[//^-dependence of Amf obtained from T = MFT over- 
estimates the charge gap obtained from the p{p) analysis 
for moderate to large coupling. The MFT is also quite 
misleading about the finite temperature antiferromag- 
netic (AF) long range order, that is necessarily absent 
in 2D. 

Spin correlations: To study magnetic correlations we 
look at the uniform susceptibility xC^) for various U/t 
shown in Fig. [U^b). We find a broad peak in x{T) at a 
temperature Tgpin below which short range AF spin cor- 
relations are important. We have also analyzed the short 
range spin-spin correlation function {St ■ Si+s) as a func- 
tion of T for various U/t and found that it too shows a 
growth in correlations at the temperature scale Tspin- 

To get a better feel for Tspi„ we also plot in Fig. [H^b) 
the uniform susceptibility for the 2D nearest-neighbor 
5=1/2 Heisenberg AF [15], that describes the low en- 
ergy spin physics of the Hubbard model for U/t ^ 1. 
Here we see a peak in x at Tspm — Jaf- To plot this 
data together with our Hubbard model results, we have 
chosen (somewhat arbitrarily) U/t = 12, so that the AF 
superexchange Jaf ~ it^/U ~ t/3. 

Except in weak coupfing U/t <g; 1, where the charge 
and spin scales are essentially identical, we find in 
Fig. [TJa) that the two scales are quite different for 
U/t > 1. In the large U (Mott) limit we expect Ach = 
U/2 > Tspin ^ Jaf = 4,t'^/U. For any U/t we expect 
to see local moments below the charge gap T < Ach, 
and the build up of AF spin correlations between the 
moments for T < Tspin. 

Entropy: In cold atom experiments the entropy S can 
be monitored more easily than the temperature. We 
calculate S{T) from QMC in two different ways. The 
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FIG. 2: (a) Density p as a function of chemical potential fj, for U = 6t and various T. Note that p is pinned at unity for 
< Ach, the charge gap. The large error bars away from half-filling are a result of the sign problem. Spin susceptibility x{T) 
for (b) positive U at half filling and (c) negative U at density p — 0.7. The peak in x{T) is used to determine characteristic 
temperature scales for the two models described in the text, (d) The superfiuid density ps{T) for various U < at p = 0.7. 
This is used to estimate the Kosterlitz-Thouless Tc in 2D using psiTc) = {2/n)T. 
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FIG. 3: Entropy per site s — S/Nks as a function oiU at a 
temperature corresponding to the characteristic scale for that 
U . (a) For (7 > 0, the entropy at Ach, the charge scale, and 
Tspin, the spin scale is plotted, (b) For U < the entropy at 
the pairing scale T = T* and the superfiuid T = Tc is plotted. 

first method [l^ II3l is to integrate down from infinite 
temperatures: S{(3) ^ pN\nA + f3E{l3) - J^dP'E{^'), 
where /? = 1/T and E is the energy. The second method 
[lit is to integrate up from T = 0. Here we fit E(T) 
data to a suitable functional form, find the specific heat 
C(T) = dE{T)/dT and calculate S{T) = J^C{T')/T'. 
The methods agree to within a few percent. 

The curves of constant S are plotted in the {U, r)-plane 
in Fig.[TJa). At weak coupling the charge and spin scales 
arc both exponentially small and the system looks like a 
highly degenerate normal Fermi gas with a low entropy. 
It is only at very low entropy (per particle) s ^ ln4 ~ 
1.386, that one observes either the charge gap or build up 
of spin correlations for U/t < 1. For moderate coupling, 
say, J7/t ~ 6 we see that s = In 2 ~ 0.693 already puts us 
below the charge gap scale. By lowering the entropy to 
s ~ 0.4 it is already possible to cross Tspin- The difference 
between the entropy required to access the charge and 
spin scales grows as U increases. We note that Tgpin sets 
the maximum scale for the AF TnccI if we were to couple 
2D layers to cut off the fluctuations and stabilize AF 
order in layered system like the parent insulator of the 
high Tc superconductor. In Fig. [31 we summarize our 
results for the entropy per site S/Nks at the charge and 
spin temperature scales plotted as function of U . 




FIG. 4: Double occupancy ini-^riii) as a function of temper- 
ature T for (a) [/ > and (b) i7 < 0. The legend in (b) 
also apply to the positive U values in (a). Also shown are 
the comparisons with dynamical mean field theory DMFT in 
panel (a) that shows a much more pronounced anomalous re- 
gion with dD/dT < 0. The inset in (a) shows the lack of 
system size dependence in the D{T) results. 

Cooling and Thermometry: We next turn to a beauti- 
ful suggestion of Werner et. al. [19] for adiabatic cool- 
ing in an optical lattice, that exploits the anomalous T- 
dependence of the double occupancy D = {rii^nii) ob- 
served in dynamical mean field theory (DMFT). In brief, 
their idea is as follows. The entropy and D are related 
to the free energy F via S = —dF/dT and D = dF/dU. 
Using a Maxwell relation we find dS/dU = —dD/dT. 
For C/ > 0, the "natural" expectation is dD/dT > so 
that the temperature increases along a constant S curve 
upon increasing U/t. Thus the DMFT observation of 
a significant anomalous region with dD/dT < implies 
that we can follow a curve of constant S and cool the 
system as the lattice is turned on, i.e., U/t is increased. 

We show in Fig. Ufa) the DMFT curve for D{T) at 
U — ()t contrasted with 2D QMC results. The absence 
of a significant anomalous regime with dD/dT < in 2D 
is likely a result of short range spin correlations that are 
important in low dimensional systems but neglected in 
DMFT. We see that the constant-S* curves in Fig. [T](a) 
do not show a significant negative slope and thus one 
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cannot obtain adiabatic cooling in 2D. 

From Fig.lDJa) we see that the double occupancy D{T) 
is quite weakly T-dependent for T < t, but its monotonic 
T-dependence at higher temperature suggests that D(T), 
that can be measured in cold atom experiments [s*], can 
be used for thermometry in the range 1 < T/t < 10. 
Such a thermometer would need QMC results for D(T) 
callibration. We show in the insert to Fig. IDJa) that the 
double-occupancy, that is a local observable, has no sig- 
nificant system size dependence and can indeed be very 
accurately determined by QMC. 

Attractive Hubbard model at arbitrary filling: 
The U > model away from p — 1 has a fermion sign 
problem in QMC, so we turn to a different model to dis- 
cuss superfluidity in a 2D lattice system. We change the 
sign of the interaction and examine the attractive Hub- 
bard model that exhibits s-wave pairing [13] and a BCS 
to BEC crossover [13[ as a function of \U\/t. We will 
consider the case of p = 0.7 filling for concreteness; we 
choose to work away from half-filling, because at p = 1 
the results would be identical to those discussed above 
using a well-known particle-hole transformation. 
Pairing (Pseudogap) scale: The analog of the charge 
scale below which moments form in the repulsive model, 
is the pairing scale T* in the attractive Hubbard model 
below which double occupancy grows. Since pairing rep- 
resents the onset of singlet correlations, one of the sim- 
plest ways of probing this scale is through the susceptibil- 
ity xC^) plotted in Fig. ^c) for a sequence of attractive 
couplings \ U\/t. The peak in x is a measure of the pairing 
scale, for below it the spin response is strongly suppressed 
by the formation of pairs. The |f7|/t-dependcnce of T* 
as determined from x is plotted in Fig. [Hb). T* essen- 
tially has the same |C/|-dependence as the T — pairing 
gap: exponentially small in \U\/t for the weak coupling 
BCS limit and proportional to U in the strong coupling 
Bose limit. For T <T* there is a pseudogap jl3i] due to 
pairing correlations in the density of states, even above 
Tc where the system is not superfluid. 
Berezinskii-Kosterlitz-Thouless Transition : At lower 
temperatures there is a transition to a superfluid phase 
in the attractive Hubbard model for p ^ 1 . There is alge- 
braic order in the pair-pair correlation function [l^ below 
Tc and a non-zero superfluid density ps . We calculate ps 
from the transverse current-current correlation function 
[20^ and use the universal jump in pa[T~) = 2Tc/7r to es- 
timate the BKT Tc as shown in Fig. ^d). At a filling of 
p = 0.7 there is a maximum in Tc ~ 0.2t as a function of 
attraction at ~ 5 as seen from Fig. [2ld). The non- 

monotonic dependence of Tc on \ U\lt is expected with an 
exponentially small Tc in the weak coupling BCS limit 
and Tc ^ /\U\ in the strong coupling BEC limit. 
Entropy, Cooling and Thermometry: We next calculate 
the entropy using the methodology described above and 
plot the curves of constant S in the (|J7|, T)-plane in 
Fig. mb). We see from Figs. [T]Jb) and 3(b) that for a 



coupling |?7|/t ~ 5 at which the Tc peaks for p = 0.7, one 
only needs to lower the entropy S/Nks < 0.8 to enter 
the pseudogap regime below T* . To actually observe su- 
perfluidity one needs to cool below the BKT transition, 
that corresponds to S/Nks < 0.1. 

The double occupancy D{T) for [/ < in Fig. IHb) 
shows a rather small regime of anomalous behavior, 
which now corresponds to dD/dT > 0. Thus the 
prospects of using adiabatic cooling [l^ in 2D do not 
look promising for the attractive case either. However, 
Fig. lUJb) does suggest that D{T) can be used as an 
effective thermometer for the high temperature range 
t<T< lot. 

Conclusions: Current fermion optical lattice experi- 
ments [4] have achieved an entropy per particle ~ In 2, 
sufficient to observe the charge gap in the repulsive Hub- 
bard model or the pairing pseudogap in the attractive 
case. Observing antiferromagnetic correlations or super- 
fluidity in 2D systems will require a further reduction in 
the entropy by a factor of three or more. It is possible 
that the inhomogeneous density in a trap can lead to 
a redistribution of entropy with some regions having a 
much lower entropy than others. 
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